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Abstract 

There exist two types of stationary instability of the flow over a rotating disc 
corresponding to the upper-branch, inviscid mode and the lower-branch mode, which has a 
triple deck structure, of the neutral stability curve. A theoretical investigation of the linear 
problem and an account of the weakly nonlinear properties of the lower branch-modes have 
been undertaken by Hall (1986) and MacKerrell (1987) respectively. Motivated by recent 
reports of experimental sightings of the lower-branch mode and an examination of the 
role of suction on the linear stability properties of the flow here we investigate the effects 
of suction on the nonlinear disturbance described by MacKerrell (1987). The additional 
analysis required in order to incorporate suction is relatively straightforward and enables 
us to derive an amplitude equation which describes the evolution of the mode. For each 
value of the suction a threshold value of the disturbance amplitude is obtained; modes of 
size greater than this threshold grow without limit as they develop away from the point of 
neutral stability. 


1 This research was supported in part by NASA Contract No. NAS1-18605 while the authors were in residence 

at ICASE, NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


The stability of rotating disc flows has been the subject of many studies, one reason for 
this being that the basic boundary layer flow given by Von Karman’s (1921) exact solution 
is fully three-dimensional. Interest also stems from the observed similarities in the type of 
instability occurring in rotating disc flow work (Gregory, Stuart & Walker 1955) and near 
the leading edge of a swept back wing, Poll (1978). This instability, known as cross-flow 
instability, is due to the inflectional character of the basic velocity profile and was first 
studied extensively by Gregory et al. (1955). They showed with china clay techniques that 
this instability took the form of a regularly spaced pattern of equi— angular spiral vortices 
which was stationary relative to the disc. Stuart (in Gregory et al . ), using inviscid theory, 
suggested that the instability could be associated with a particular inflectional profile in 
which the inflection point coincided with a point of zero velocity somewhere in the flow. His 
calculation gave the predicted number of vortices to be approximately four times greater 
than the observed value of about 30 but the angle of 13° between the axes of the vortices 
and the radius vector was in excellent agreement with their experiments. 

Several attempts have been made to explain the observed spiral patterns by means 
of a linear stability theory (see, for example, Cebeci & Stewartson (1980), Malik et al. 
(1981), Mack (1985), Malik (1986)). Malik et al. (1981) showed that Coriolis force and 
streamline curvature cannot be neglected- a point first made by Faller & Kaylor (1966) and 
Lilly (1966) in the context of the Ekman boundary layer. The theoretical results of Mack 
(1985) strongly suggest that the spiral streaks observed in flow visualisation experiments 
are the constant phase lines of the merged wave patterns produced by several random 
sources on the disc. This was first found by the experimental investigations of Wilkinson 
& Malik (1985) who revealed that the wave patterns from each point source spread out 
circumferentially downstream of the source and that the wave patterns from a number of 
sources eventually merge and cover the entire circumference of the disc. 

Experiments by Fedorov et al. (1976) showed that for various rotation rates of the disc 
an instability occurred which also appeared as a pattern of spiral vortices which numbered 
in the region 14 to 16 and which had their axes inclined at angles of about 20° to the 
radius vector. These discrepancies between these observations and those of Gregory et al. 
(1955) imply that the vortices seen by Fedorov et al. were not those of Gregory et al. 

Malik (1986) calculated the neutral stability curve for stationary disturbances and 
demonstrated that two types of mode could exist. The first, the upper branch of the 
neutral curve, corresponds to the inviscid mode described by Stuart in Gregory et al. The 
second, the lower branch mode, is an essentially viscous disturbance which corresponds to 
zero wall shear stress of the crossflow velocity profile. 
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A linear asymptotic investigation of the inviscid mode at high Reynolds number was 
conducted by Hall (1986); work that has been extended to include nonlinear effects by 
Gajjar (personal communication). Further, Hall (1986) demonstrated that the lower 
branch viscous mode is governed by a triple deck structure (which is derived from the 
classical structure as reviewed by, for example, Smith 1982) and it is this mode which 
is thought to have been observed by Fedorov et al. (1976). MacKerrell (1987), hereafter 
referred to as M, extended the study of Hall (1986) in order to give an account of the weakly 
nonlinear properties of the lower branch mode. Following the framework set up in Hall 
(1986) M derived solutions which enabled the finite amplitude growth of a disturbance close 
to neutral to be described. Importantly, she found that nonlinearity has a destabilising 
effect on the vortices and discovered the existence of a threshold amplitude (in the sense 
of Stuart 1960, 1971 and Watson 1960) such that disturbances smaller than the threshold 
value decay as they move away from the neutral position whereas those above threshold 
grow explosively. This led M to conclude that for small disturbances the inviscid mode 
of Gregory et al. (1955) is probably dominant whereas for larger amplitude perturbations 
the short wavelength mode observed by Fedorov et al. (1976) is the more important. 

Until very recently there have been, to the best of the authors’ knowledge, few, if 
any, other observations of the short wavelength disturbance. At the 1991 ICASE/LaRC 
workshop on transition and turbulence Corke reported that in his experiments on rotating 
disc flows a subcritical stationary instability had been observed- the mode which is almost 
certainly that seen by Fedorov et al. (1976) and analysed in M. At the same meeting 
we learned of an investigation by Dhanak (1991) who has extended the calculations of 
Malik (1986) to include the effects of suction. The role of suction in practical flows is an 
important one for it is believed that imposition of suitable suction on a boundary layer 
can often stabilise the flow and thence delay the onset of transition. Dhanak (1991) found 
that if suction is applied to the rotating disc flow then the overall effect, at least according 
to linear theory, is that the flow is indeed stabilised: conversely blowing destabilises the 
flow. This finding invites the question as to the effect of suction on the nature of the short 
wavelength mode described by Hall (1986) and M and it is this that we address here. 

In §2 we follow the nonlinear analysis of M and indicate the modifications required in 
order to incorporate suction into the flow. We then execute some numerical work to obtain 
a quantitative description of the flow properties and conclude with some discussion. 
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2. Formulation of the problem 

We consider the case in which the disc z = 0 rotates about the 2 -axis with angular 
velocity ft. Relative to cylindrical polar co-ordinates (r, 6 , 2 ) which rotate with the disc, 
the continuity and Navier— Stokes equations, when suitably non— dimensionalised, become 


V.u = 0, 


(2.1) 

<9u 

+ (u.V)u + 2(k A u) — rr = - 

Vp+^u. 

(2.2) 

where u = ( u,v , w) are the velocity components, p is the fluid pressure, r and k are unit 
vectors in the r— and 2 — co-ordinate directions and the Reynolds number Re = Q.L 2 /u 
with v the kinematic viscosity of the fluid and L a reference lengthscale. The Reynolds 
number is taken to be large throughout. 

For the ensuing analysis it is convenient to define e = Re~& and with 2 = e s r] the 
basic steady flow is given by the solution 

u = ub = {ru(ri),rv(ri),e 8 w(r})') , 

P = PB = p{rj), 


where u, v, w and p satisfy 



u 2 — (1 + v ) 2 + u'w = u" , 2u(l + v) + v'w = v" , 

(2.3a, b) 

w' + 2u = 0, p' + ww' — 

w" = 0. 

(2.3c, d) 

The appropriate far field conditions take the forms 



u — > 0, v — ► —1 as 

v — ► 00 , 

(2.4a) 

whereas on the disc surface we impose 



u = v = 0, 


(2.46) 

together with 

w = —5, 


(2.4c) 


where S is a measure of the suction ( S > 0) or blowing (5 < 0) applied to the flow. 

The linear stability of the basic flow (2.3), (2.4) in the absence of suction or blowing 
was addressed by Hall (1986). He demonstrated that the flow is susceptible to a stationary 
disturbance which is governed by an appropriate triple deck structure whose lower, main 
and upper decks are of thicknesses 0(e 9 ), 0(e 8 ) and 0(e 4 ) respectively. Hall showed that 
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the disturbance wavenumbers in the r and 6 directions are 0(e 4 ) and solutions were 
sought proportional to 


E = exp 





(2.5) 


where the wavenumbers a and 3 expand as 


a = qq + e 2 c*i + e 3 a-i + . . . , (3 = /?o + e 2 /?i + e 3 /?2 + (2-6) 


It was found by Hall (1986) that, as shown by Gregory et al. (1955), the ‘effective 
velocity profile for a three-dimensional disturbance with wavenumbers a and (3 in the r 
and 0 directions is rau + (3v. However, if the effective wall shear rau' + (3v' does not vanish 
the modes are necessarily time dependent for a and {3 real. Therefore, in order to seek 
stationary modes which are neutrally stable at the location r the values atj and j3j must 
be chosen so as to reduce the effective wall shear to zero at leading order. 

MacKerrell (1987) extended the work of Hall (1986) to a weakly nonlinear setting. She 
demonstrated that in order to obtain a classical evolution equation for the disturbance it 
is necessary to examine disturbances of amplitude 0(e%) within the lower deck. Further, 
the mode then evolves in an 0(e) neighbourhood of the neutral position r = f and, as in 
usual nonlinear triple deck work, it is the lower tier of the triple deck where nonlinearity 
manifests itself. Our objective here is to generalise the work of M to investigate the role 
of the suction parameter S and in the following we deliberately give the briefest of details. 
The analysis differs only slightly from that of M and so we highlight the main differences 
between the present work and that contained in M. The omitted details may be found in 
M and the reader is referred to that paper where necessary. 

Following the argument of M we define a new variable r\ in the neighbourhood of the 

neutral point f by 

r = r + eri, (2.7a) 


and seek disturbances whose radial dependence may be written in the form 


d_ 

dr 


ido la i 

T~ 


i 

+ - 


e dr\ 


(2.76) 


Within the lower deck, where we define the 0(1) co-ordinate f by 2 = £ 9 £, we suppose 
that the basic flow velocity components (2.3) expand according to 


u =euo£ + £ 2 ^i^ 2 + e Z n 2 e + • • • ? 

v =evo£ + + e 3 u 2 £ 3 + • • • > (2-8) 

w = — S + e 2 uh£ 2 + e 3 W2i 3 + 
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If we perturb the basic flow by writing 


( u , v , w,p) = (u + U,v + V,€ S w + W,p + P ), 

then in the lower deck the disturbance quantities develop according to 
'rj 0C1 


(2.9) 


V=ei 


(uq + 2eU\( ) + . . .) 4 b U 0 + . . . 


E 


L £ 0o 

+ harmonic terms + complex conjugate, 
together with a similar expansion for V. In addition 

ijo£ 5 C\ 


(2.10a) 


W =£ = 


00 


[(a 0 ru! + 0oui)f 2 + •••] +e 6 Wi + e 7 W 2 + 


E 


(2.106) 


+ harmonic terms + complex conjugate, 


and 


P=e 7 > [{e z P x +e A P 2 + ...)] E 

+ harmonic terms + complex conjugate. 

Here C\{r\) is a scaled amplitude of the disturbance and we have made the definition 


(2.10c) 


-2— 2 , 00 
7o = a 0 + 


( 2 - 11 ) 


These forms of the disturbance, in which all the unknowns are functions of £ and r x , were 
first derived by M who also elucidated the forms of the harmonic terms. However, we find 
that the crucial harmonic and mean flow terms which interact to drive the fundamental 
and hence lead to the desired amplitude equation are unchanged by the imposition of 
suction or blowing. Hence, for brevity we do not repeat these results from M and now 
merely state how her overall findings need to be modified. 

Firstly, following the procedure adopted in M, we find that W 2 satisfies 


ds 3 


2 - s + 2sW 2 =A"*7' ]D x - 2r 1 7 0 2 C 1 ^- +i7oSo0iA-5 ' s2riCl 


ds 2 


0o 2 r 
_i_ C x s 2 


+ iki(raiu 0 + f3iv 0 )A *s 2 -i(3 2 v 0 A 4 7 o 


00 


_ 2 _ ~*A*s 2 _ uqA 

— l*c 2 c x 2i7 0 rxCi 


0o 4 


rK) 


— ruoA 4 


0o r 

_i 7 o s 2 dC x * i loD\ 


v 2 \ U(Q,V2s) 
U( 0,0) 


0o dr i 


+ 2SA* 


0o 


+ RHS, 


(2.12a) 
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where A = i(aofui + flovi), Di(ri ) and C 2 (ri) are functions of C\ given in M, A7 is a real 
constant whose precise value will turn out to be unimportant, s = A*£ and an asterisk on 
a quantity denotes the complex conjugate of that quantity. Further, U(a,x) denotes the 
usual parabolic cylinder function (see, for example, Abramowitz & Stegun (1964)) and 


_ 2 

RHS = — 2ir^-Uo&-% — ikiA~*s 4 (raoU2 + P 0 V 2 ) + 2*7oCi 


A' 1 _ _ U(0,y/2s) 


u — — j, — w ^ . — u ~ x 1 1 — ^ 

i(raiu 0 + Piv 0 )A~% + 3iA~ 2 s 2 (ra 0 U2 + /?o ^2 ) j 


-SUqVo 


To C.F.W + ^C, 


+ 


(l + t 

V «0 

i(ra 0 u 2 + j3 0 V2)A~ 2 s z + i{ra\UQ + ^xv 0 )A~^s 

^ w + f Cl ( 1 + |) 


\ u 0 F 2 (s) 

J U( 0,0). 


upF^if) 

U(0,0) 


S A 


-1 


tfCiFlXs) + 


(-1) 


uoFZ(s) 


U( 0 , 0 ) 


where F\(s) and F 2 (s) satisfy 

F{" - s 2 F[ + 2 sFj = 1, 


F 1 (0) = F 1 (oo) = 0, 


^( 0 , 0 ) 


(2.126) 


and 

F"' - s 2 F! i + 2 sF 2 = F(0, V 2), F 2 (0) = F 2 (oo) = 0. 


The solution for Uq when 5 = 0, say Uqo, is given in Hall (1986). It is then straightforward 
to obtain 


Uo = Uoo+SA~*C 1 r 


7qU 0 sU(0,s\/2 ) 

V2P 2 U( 0,0) 


In order for the disturbance to vanish on the disc surface we find from the continuity 
equation that dW 2 /ds = 0 on s — 0. By solving (2.12) subject to this condition we obtain 
the desired amplitude equation in the form 


dC 1 
dr 1 


(a + ic)r, C, +(!, + id)C, |C, | 2 + (e + i/)C, , 


(2.13) 


where a, 6, c, d, e and / are real valued constants, the more important of which are detailed 
presently. Multiplying (2.13) through by C*, adding the complex conjugate and replacing 
77 by r[ = 77 + (e/a) yields the amplitude evolution equation 


^■(|C 1 | 2 ) = 2or' 1 |C,| 2 +26|C 1 |*, 


(2.14) 
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and the form of the solution of (2.14) depends on the signs of a and b. ^From (2.3) and 
(2.4) we find that 

1 Suq . _ Sv o 

u i = and V\ = — 

1 2 2 2 

and recalling the imposition that the leading order wall shear (= Oi^fUQ 4- Po ! -'o ) must be 
necessarily zero in order to obtain stationary vortices we have A = —if a o/2. Following 
the manipulations of M we finally obtain 

Plloh 

a = — — -r-j 

r 2 u 0 (aor)? 



and 


where 


_ 3(o 0 f) 2 7q v$\ 

b ~ 2V2 rft V + *l)»' 


(2.15 b) 


A 2 


and 


* 

470)80*1 ~ ^ ( fa i“o + 0 i vo ) 


+ 


_ 1 


7o0o*i - (-y-) ? (raiUo + 0i « 0 ) 




(2.15c) 

h = f?( 0) = 

f 6Ui f' d \dd = 0.59907, 
h 217(0,0) 

(2.15d) 

, _ 2*5(0) 

r°vy,0) de = 0A5m5 

(2.15e) 

17(0,0) 

Jo ^ 2 ( 0,0) 


Consequently, we now have the outcome that the formulae (2.15) for the coefficients of 
the amplitude equation appear to be independent of the suction parameter S and hence 
are identical to those of M. However, the constants appearing in the expansions (2.8) for 
the basic flow quantities close to £ = 0 are functions of S and thence so are a and b. In 
the following section we conduct a few calculations for various values of S and draw some 
conclusions. 


3. Results and discussion 


In order to calculate the coefficients a and b in (2.14) we merely need to solve the 
basic equations (2.3) subject to boundary conditions (2.4) and note that the function 
(fai^o + 0i vq) appearing in definition (2.15c) is given by 


7o*i + 


vy Q u Q 

“0T 


1 + rf ) *2 = *A* (ouruo + 0i wo) -Q, 

Po 
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a relation which was deduced in M. This relation also serves to determine the value of 70 and 
this, in addition to the requirement that the leading order wall shear aruo + (3 oVq vanishes, 
furnishes the leading order radial and azimuthal wavenumbers ao, /?o- The equations for 
the basic flow were solved using a suitable NAG routine and then we found it convenient 
to write 


a = Ar * , b — Br ^ , 


(3.1) 


where the functions A(S) and B(S) are given in Figure ( 1 ). We note that for all values of 
suction or blowing, — 00 < S < 00 , we have A < 0 and B > 0 and the former condition 
implies that the amplitude of the solution increases or decreases depending on whether r 
is less than or greater than the neutral value. It is also straightforward to demonstrate 
that if we write 

|Ci | 2 = (^y/—aj2b^ y and rj = x/y/^a 

then (2.14) becomes 


j- = -2 xy + y 2 , 
ax 


whose solution is 


y(x) = 


(» " *r eTf ( x )) ’ 


where yo = y(0) and erf(x) = JJ e -t dt. Since y(x) must be greater than or equal to 
zero and 0 < erf(x) < 1, the solution is bounded and valid for all x if yo < 2/y/n. However, 
for yo > 2/v/SF, y(x) becomes infinite at some point xo, 0 < xo < 00 , and the solution is 
only valid for x < xo- Thus the value yo = 2j y/n represents a threshold between solutions 
which grow indefinitely and those that decay to zero. The corresponding disturbance 

threshold value 

. 2 I cT 2 / A _ 3 

Cle= V*v 2& = (3,2) 

is shown in Figure (2). We note that for large suction ( S 1) the threshold amplitude 

becomes small whereas for large blowing this amplitude grows. We now verify these 
qualitative results by considering the appropriate limits. 


3.1. The large suction limit, S 1. 

For 5 » 1, Stuart (1954) has demonstrated that the basic flow profile is confined 
within a zone of thickness 0(1/S) and if we write £ = 5 ? 7 where £ = 0(1) then, to leading 
order, 


u = 


2 S 2 


(< 


-< - 


2< *) , v = 1 — e ^ 


and w = —S. 
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Consequently the constants defined in the basic flow expansions (2.8) assume the forms 

1 - c 

U ° = 25’ V ° = ~ S 

and these lead to the asymptotic expressions 

A 2%I$I~*S¥ « -1.20785^, (3.3a) 

i 

24 7 4 -77 - 77 

B i r 5 T « 8.19525^, (3.36) 

(21 - 10 V 2 )/ 2 7 

where /i and I 2 are given by (2.136, c) and A and B are the coefficients in the amplitude 
equations (2.14) and which were defined in (3.1). In turn, the threshold amplitude 

|Ci| e « 0.30635"M. (3.3 c) 


3.2. The large blowing limit, S ■< — 1. 

In his consideration of boundary layer flows Pretsch (1944) demonstrated that the 
influence of viscosity becomes almost negligible for sufficiently large blowing and this work, 
together with the results emanating from studies of fluid injection in supersonic boundary 
layers by Smith Sz Stewartson (1973a, b), enable us to derive the flow structure for the 
present problem when S <C — 1. We find it convenient to define the positive quantity 
S = —S. It is now the case that the extent of the boundary layer increases and formally 
becomes O(S). Therefore, if we define the 0(1) co-ordinate z\ by z\ = rj/S we find that 
in the region 0 < Z\ < it / 2 the basic flow solution develops according to 

u = Fq(zi)+^F\(zi) + . . . , v = Go(^i) + -^jGi(zi) + • . . , 

S . S (3-4) 

i» = Sffo(zi)+ jffi(zi) + ..., 

with 

Hq = Go = cos 2 z\ and Fq = coszi sin^j. (3-5) 


The profile decays algebraically as z\ — > 7r/2 and a thin viscous layer centred on this 
point ensures that the solutions for u and v are brought to zero at infinity. Details of this 
viscous structure are not needed for present purposes as we are solely concerned with the 
values of u'(0) and u'(0). Consideration of the problem for G\(z \ ) enables us to conclude 
that the flow constants uq and uq in (2.8a, b) are given, at leading orders, by 



and 
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In turn this leads to the asymptotic results 


A ~ -2*I 2 * Jj 4 S'* w -0.3020 (-5) 4 , (3.6a) 

I 

6J 4 n. _ _ ii 

5 ! r 5 T » 2.0488 (-5) 4 , (3.66) 

(21 - 10\/2)J 2 5 

with a threshold amplitude given by 

\Ci\ e « 0.3063 (-5) ^ ft. (3.6c) 


4. Conclusions 

On Figures (1) & (2) we have indicated the asymptotic behaviours (3.3) & (3.6) both 
for 5 < —1 and for S > 1. We observe that these one term asymptotic predictions are 
very accurate for a surprisingly large range of values of S ; indeed they are graphically 
indistinguishable from the respective curves for |<S| greater than about five. 

The most important deduction to be drawn from Figure (16) is that the nonlinear 
coefficient in evolution equation (2.14) is positive for all values of blowing and suction. 
Thus in all cases nonlinearity has a destabilising influence and we have found the existence 
of a threshold amplitude |Ci| e as discussed by Stuart (1960, 1971) and Watson (1960) such 
that all disturbances of amplitude less than threshold decay whereas those greater than 
threshold become infinite as the distance from the neutral stability position increases. This 
in turn leads to turbulence. 

Our results described above seem to suggest, perhaps unexpectedly, that the effect of 
including suction in the rotating disc problem tends to lower the threshold amplitude for 
the disturbance, see Figure (2) and result (3.3c). Moreover, Figure (2) indicates that only 
very moderate suction or blowing is needed in order to have a significant effect on the 
critical disturbance amplitude. In the work of Dhanak (1991) it was found that suction 
stabilises the flow in as much that inclusion of suction reduces the region of wavenumber/ 
Reynolds number parameter space in which the basic profile is unstable to stationary, 
infinitesimal modes. However, since the threshold amplitude decreases with increasing 
suction we speculate that an experimental configuration would need to be less strongly 
forced in order for the subcritical instability to occur than would be the case for the zero 
suction problem:- conversely, imposition of blowing destabilises the flow according to linear 
analysis but the much larger threshold amplitude suggests that exciting the subcritical 
disturbance experimentally might be extremely difficult. 
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We now say a few words about the extension of our present problem into the 
compressible regime. Seddougui (1990) considered the effects of compressibility on the 
subcritical mode of M and found that for both an adiabatic wall and an isothermal wall 
the stationary mode is only possible over a finite range of Mach numbers. The analysis 
conducted by Seddougui (1990) is naturally more involved than that of M but we have seen 
here that, at least for the incompressible flow problem, the inclusion of suction requires 
only relatively minor changes to the nonlinear analysis. With suction the variations in 
the coefficients of amplitude equation (2.14) are predominantly due to the changes in the 
wall shears; these are in turn due to the changing boundary condition (2.4c) for the basic 
flow. Preliminary consideration of the compressible version of the present problem also 
suggests that the modifications required in order to account for suction in Seddougui’s 
(1990) work will also be relatively straightforward. Further studies on this aspect are in 
progress together with an investigation into the properties of the time dependent version 
of the mode considered here. 

Finally there is the question of the extension of the present work into the fully nonlinear 
regime. Examples of fully nonlinear calculations for suction/blowing problems are typified 
by those of Smith & Stewartson (1973a, b) for plate and slot injections into supersonic 
boundary layers. Extension of our work along the lines suggested in these papers would 
be valuable indeed. 
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logio(-A) 


Figure (1). The dependences of amplitude equation coefficients A and B defined by (2.14), 
(2.15) and (3.1) on the suction parameter S. a) log 10 (— A): dotted lines indicate the 
asymptotic values (3.3a) and (3.6a), b) log l0 B: dotted lines indicate the asymptotic values 
(3.36) and (3.66). 
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